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Abstract 

We describe a way to continue the fermionic renormalization group flow into 
phases with broken global symmetry. The method does not require a Hubbard- 
Stratonovich decoupling of the interaction. Instead an infinitesimally small symme- 
try-breaking component is inserted in the initial action, as an initial condition for 
the flow of the selfenergy. Its flow is driven by the interaction and at low scales 
it saturates at a nonzero value if there is a tendency for spontaneous symmetry 
breaking in the corresponding channel. For the reduced BCS model we show how 
a small initial gap amplitude flows to the value given by the exact solution of the 
model. We also discuss the emergence of the Goldstone boson in this approach. 

1 Introduction 

Renormalization group (RG) flows for interacting fermions have experienced a mas- 
sive upswing in recent years [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]. Besides their conceptual 
importance for the understanding of many-particle systems they have become a suc- 
cessful unbiased vehicle for the detection and evaluation of Fermi surface instabilities 
in low-dimensional many-fermion systems. For example, the two-dimensional Hub- 
bard model has been subject to numerous RG studies over the last years. For most 
parameter regions in the weak coupling range of this model, it has led to a good qual- 
itative picture of the dominant instabilities of the Fermi liquid state. In most cases 
the phases suggested by the RG are states with spontaneously broken symmetry, like 
superconducting or magnetically ordered states. Yet a major drawback of the existing 
weak-coupling fermionic RG schemes is their failure to give a controlled access to 
the symmetry-broken regime. In the language of the fermionic RG used up to now, 
the instabilities towards these states give rise to a flow to strong coupling at some 
small energy scale where at least one class of coupling constants grows larger than 
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the bandwidth of the system. Then the weak-coupling RG becomes unreliable, and if 
one ignores this and tries to continue the flow, the coupling functions diverge before 
all degrees of freedom have been integrated over. Hence, so far the only way to de- 
velop a theory for the strongly coupled phase is to stop the fermionic RG flow before 
it becomes unreliable and to resort to some other method for the scales below, e.g., a 
mean-field treatment, bosonization [12, 13, 14], or exact diagonalization of a heavily 
reduced model [15]. The antiferromagnetic phase of the half-filled Hubbard model 
was studied with bosonic RG methods in [16]. 

Here we propose an extension of the fermionic functional RG methods which al- 
lows us to access the symmetry-broken regime continuously. The basic idea is to insert 
a small symmetry-breaking component into the initial condition of the flow at the ini- 
tial energy scale eo- ln a magnetic system, this would be a small external field, and if 
spontaneous symmetry breaking happens, the spontaneous magnetization is obtained 
in the limit where this field vanishes. For the case of a Cooper instability towards a 
superconducting state the symmetry-breaking term is a small superconducting gap am- 
plitude Ao. In the RG flow, degrees of freedom with energy above a scale e s = eoe~ s 
are gradually integrated out, and the full model is recovered in the limit s — > oo, where 
e s — > (eo is a fixed energy scale; we take it equal to the bandwidth of the model). Now 
in the course of the flow, the gap amplitude gets renormalized by the flow equation for 
the self-energy, so that A s tends to a (A -dependent) final value A nnal in the limit of 
vanishing energy scale, s — > oo. We will show below that, together with a reorganiza- 
tion of the hierarchy of flow equations for the one-particle irreducible vertex function 
proposed recently by Katanin [17], this scheme allows us to reproduce the exact mean- 
field results for the reduced BCS model. That is, we can take the limit s — > oo, hence 
integrate over all degrees of freedom without encountering a divergence in the cou- 
pling functions, and when the external field Ao is sent to zero after that, the final gap 
Afi na i converges to the solution of the BCS equation. Moreover, the two-particle ver- 
tex also converges to the exact result for the BCS model. Furthermore we show that 
the Ward identity from global U(l) charge symmetry is satisfied in our approach. This 
Ward identity, together with the solution for the gap, fixes the two-particle vertex of 
the model, and we recover the divergence at zero energy which is the remnant of the 
Goldstone mode in the reduced BCS model. 

In this paper we restrict our detailed calculations to the reduced BCS model, but we 
emphasize that our method is also applicable and expected to give reasonable results 
for models where no exact results are available and where no mean-field theory can be 
developed straightforwardly. This will be discussed in Sec. 5. 

The structure of the paper is follows. In Sec. 2 we briefly go through the essential 
points in the RG formalism for the one-particle irreducible (1PI) vertex functions the 
further treatment is based on. We describe a modification of the scheme proposed 
recently by Katanin[17] which allows one to reproduce the selfconsistent random phase 
approximation (RPA) + Hartree resummation in a certain channel. In Sees. 3 and 4 we 
apply this modified scheme to the reduced BCS model and show how it allows us 
to flow into the symmetry-broken regime. We also show how the Ward identity is 
recovered for this case. The outlook in Sec. 5 briefly discusses the superconducting 
gap flow in models more general than the reduced BCS model. The appendix contains 
some details on the derivation of the gap flow in Sec. 4. 
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2 The renormalization group flow 



2.1 The flow for the 1PI functions 

We briefly recall the flow equations for the 1PI vertices of a general fermionic model, 
as described in [4] (for setup and notations, see [4]). The generating functional W(H) 
of the connected Green functions generated by sources H is defined as 



e 



-W(Q,H) = J c -2(*>Q*)-Vo(*) + (ff,*), (1) 



where the quadratic part of the action, Q, defines the free propagator and we assume 
that the interaction Vq is short range and even, Vq{— v &) = Vo(^>). We have explic- 
itly denoted the dependence of W on Q here. Most of the time we shall abbreviate 
W(Q, H) = W(H); the Q-dependence will resurface in the notation when we dis- 
cuss the Ward identities. 

In an RG flow, Q = Q s depends on a scale parameter s. Hence W also depends on 
s and it satisfies the differential equation 1 (see [4]) 

= -'I (si '« mY w (2 > 

Here the dot denotes the derivative with respect to s. The Legendre transform T of W 
also depends on s, and it satisfies the differential equation (see [4]) 



f(*) = i(tt,Qtt) + ltt 



(3) 



This nonpolynomial equation can be transformed into a hierarchy of differential equa- 
tions in which the vertex functions appear only polynomially by expanding in the fields: 

r(#) = iv + i(*,f (2) *) + r^ 4) (i') (4) 

where 

r ( - 4) (*) = J2 r (2m) (*) (5) 

m>2 

and r^ 2 " 1 ' is homogeneous of degree 2m in "J. There are no odd terms in this expansion 
because V is even in As in [4] we denote the derivative appearing in (3) by 

Txx '- 5* X 5* X ,-- Txx ' +Txx ' (6) 
For V = 0, f( 2 ) = Q. For V ^ we define the selfenergy S by 

f (2) = Q - S. (7) 

'The form stated in (37) of [4] follows by straightforward differentiation, (up to an inessential field- 
independent term arising from the different normalization of the Grassmann Gaussian integral in the defini- 
tion of W). 
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The full propagator is 

G= (f^y 1 = (Q-S)- 1 . (8) 

With this, 

f- 1 = - (l-Gf (2;4) ) G. (9) 
Comparing homogeneous parts in we arrive at the set of equations 

K = -^Tr(QG) 

t = ~Tt{ST^) (10) 
p (2m) = I Tr (5f( 2ro + 2 )) + ip 2m Tr(sf^ 4 )Gf^ 4 )(l-Gf^ 4 ))- 1 ) 
where 

5 = -G(3G (11) 

is the single scale propagator and P2 m denotes the projection on the degree 2m part of 
any element of the Grassmann algebra. Because every factor f(- 4 ) increases the power 
of \I> by 2, only finitely many terms appear on the right hand side of the equation for 

f (2m) equation foj . y = p(4) readfj 

V = ^Ti(S£W) + ^Tt (j5VGVy (12) 

For m = 3 we have 

f(6) = I-ft(sfW) (13) 
i (Tr(5f (6) GV r ) +Tr(SVGf {6y ) +Tr(SVGVGvj) . 



+ 



In general, the flow equations have the property that on the right hand side of the 
equation for V, a term quadratic in V appears, but that for all m > 3, r( 2m ) appears at 
most linearly on the right hand side of the equation for T^ 2m ^> . The term linear in r^ 2 " 1 ) 
leads to a resummation of four-point insertions which includes ladder summations, as 
described in [4]. 

It is thus natural to count powers in an expansion in the renormalized coupling 
function V. We assume that the initial interaction at s = consists of only a selfenergy 
term and a four-point vertex Vb- By forming a tree and adding one extra loop line we 
see that a 1PI graph with r four-legged vertices can have at most 4r — 2(r — 1) — 2 = 2r 
external legs. Because I^ 2 " 1 ) is a sum of 1PI graphs with 2m external legs, it is at least 
of order m in V. Thus, to order m in V, we can drop r^ 2m+2 '. Then the equation for 
r( 2m ) closes. This procedure is described in detail for the second and third order flow 
in Section 2.2; it generalizes straightforwardly to arbitrary orders. 

For an antisymmetric operator A — —A T , define the functional Laplacian by 

A -< - (w A <k) ^ l AX S AY mx) A{x ^wm- <14) 
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Laplacians can be used to write the traces in a convenient form; we shall in the follow- 
ing use that 

Tr ( 5 f(2m+2)) = _A s r( 2m+2 >. (15) 

The minus sign appears because S is antisymmetric. Representations by Laplacians can 
also be used to rewrite the other terms in the flow equation, by introducing copies of the 
fields. This will be useful in the derivations of the next section, so we discuss it here. 
We introduce two (or more) copies , , ... of the field These copies serve to 
simplify applications of the product rule and the distinction of which derivatives act on 
which function. Since the copies are only a combinatorial device, we shall in the end 
evaluate at "J/W = = ... = *. We denote this operation by [ • ], e.g. 



^(^(i))B(^(2))c($(3))| = A(*)B(*)C(tf). 



With this and the notation 



we can write e.g. 



A W) = ( _J_ C 6 



Here we also used the abbreviation = V(\&W). 



(16) 



(17) 



(18) 



2.2 The second and third order G scheme 

In this section, we justify the modification proposed in [17], namely putting G propaga- 
tors instead of single scale propagators S on the lines in the 1PI hierarchy, and describe 
a systematic procedure for doing this. In second order, this replacement is simply 



A g, 2) A U,2) y (i V (2)j = _ 1 | (^( A g' 2) ) 2 ) V^V^ 



(19) 



We show below that this suffices to do flows with symmetry breaking that recover the 
BCS solution exactly. However, a further discussion is necessary because this replace- 
ment is to some extent ad hoc. Indeed, we shall see that its systematic justification 
requires going to all orders in the expansion in the renormalized coupling. To any fixed 
order in V, the replacement is not unambiguous because the difference G — S = GSG 
is at least of first order in V. Thus, in any contribution of order m, a single-scale 
propagator S associated to a line can be replaced by a differentiated full propagator G 
up to orders V rn+1 - or one may as well leave it as a single-scale propagator S up to 
orders v m+1 . The replacement of S by G may be more suited for self-consistency 
arguments, but in a treatment where powers of V are kept only up to a fixed order m, 
one cannot chose one over the other. The ambiguity is removed if one realizes that one 
can rewrite the 1PI hierarchy to all orders in a way where only full four-point func- 
tions V appear as vertices, and full propagators G and differentiated full propagators 
G appear on the lines, and where all tadpole terms are removed, except in the equation 
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for E itself (i.e. the lowest equation in the hierarchy). We give the explicit equations 
to third order in V. They have a simple and nice form which will be useful for doing 
flows which fully take into account the vertex corrections that one needs to do in a full 
two-loop calculation of the selfenergy flow, which has not been done up to now. 

As already discussed above, we assume in the following that the initial conditions 
for the flow are that T^- 6 ) = at s = 0, and that G(0) = because the cutoff function 
is such that s = corresponds to energies outside the band range. 

The derivative of (8) gives 



G = S + GEG. 



(20) 



The initial condition for is zero, so (13) implies that r' 6 ^ is at least of third order 
in V. The terms in which and V appear in (13) are therefore at least of fourth 
order, and to third order, the equation for simplifies to 



f( 6 ) = ^Ti(SVGVGV). 



(21) 



By (20), replacing S by G in this equation does not change anything to third order 
since £ is at least of first order in V. By (12), V is at least of second order in V, so 



1 d 



f (6) = Tr(GVGVGV). + O (V 4 ) 
6 os 



(22) 



Discarding the fourth and higher order terms we can integrate this equation. Because 

G(0) = and r< 6 ) (0)^0, 



T {6) = ^Tr(GVGVGV). 



Inserting this into (12) and using (15), we get 



V 



Tr (sTGf) 



\& s \mgvgvgv) 



(23) 



(24) 



If both derivatives of | As act on the same factor V, the result is t^AsV which equals 
S. There are three such terms, and they all give the same contribution by cyclicity of 
the trace. Thus the 1/3 gets canceled and by (20), we can combine this term with the 
second order term in (24), to replace S by G. Using (18), we obtain the flow equations 
that are exact to third order in V as 



V = -- 



lw f d 



<9s 



ld_ 

49^ 



Ag' 2) Ag' 3) (Ag< 3) ) 2 VWVWV® 



E = 



1 



-A S V. 



(25) 



(26) 



The change in third order, as compared to (19), is simply the s-derivative of a single 
third order diagram. One can also rewrite the equation for the selfenergy such that no 
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Figure 1: The graphical representation of the third order flow equations. The dot de- 
notes the derivative with respect to the scale parameter s. All lines in the graphs carry 
full propagators G, except for the one with a dot, which carries a G, and the tadpole 
line in the equation for S, which carries a single-scale propagator S. 

^-propagator appears any more: by (20), (26) becomes 

(*,±V) + A GtG V = -^A 6 V (27) 

If we can solve these equations, the 1PI six-point function is given by (23) and I^ 2 ™- 1 = 
for all m > 4. That is, the connected Green functions with 8 or more external legs 
are given by tree graphs made from V and T^. 

If one wants to count loops instead of powers of the renormalized coupling func- 
tion, this system of equations is exact to two-loop order. However, we believe it is 
more natural to think of an expansion in V than a loop expansion. 

The action of Ag can be described graphically as contracting any two of the exter- 
nal lines corresponding to a factor "Jy to an internal line, thus reducing the number of 
external legs to four. Terms where both derivatives of a Laplacian act on the same ver- 
tex factor V correspond to tadpole diagrams. If the derivatives act on different vertex 
functions, we get a line of the graph connecting these vertices. This gives the graphical 
representation of (25) and (26) depicted in Figure 1. 

An important feature of (25) is that the product of propagators G is differentiated. 
This will make it possible to integrate the flow equation when only particular quadratic 
terms are kept on the right hand side. This corresponds to graphical resummations - the 
RPA if only the particle-hole RPA-type graph is kept, and the BCS ladder summation 
if the particle-particle terms are kept. A similar structure appears in the Wick ordered 
form of the RGDE (see [18], Section 4.5.4 and 4.5.5). However, in (25) the selfenergy 
is taken into account, which was not the case in the corresponding formulas of the Wick 
ordered flow. 

2.3 Self-consistency 

In this section we illustrate how such a resummation works at the example of the RPA 
plus Hartree resummation, which was originally given by Katanin [17]. 

We drop the third order term from (25). The remaining quadratic term in V is 
identical to that in (52) of [4] except that the propagators S and G are replaced by 
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an s-derivative of two propagators G. Thus all considerations about symmetries and 
the derivation of the equations for the coefficient functions remain unchanged, and the 
equations (88)— (91) of [4], as well as their graphical representation in Figures 4 and 5, 
carry over unchanged. The only replacement is that instead of (92) we now have 

C±(q,k) = -?- s (G(k)G(q±k)). (28) 

The RPA plus Hartree resummation corresponds to keeping only the RPA type particle- 
hole graph in Figure 4 (c) and only the Hartree graph in Figure 5 of [4]. This corre- 
sponds to the equations 

V(pi,P2,P3) = J dk C+(p2-p3,k)V(k,p 2 ,P3) V(p 1 ,p 2 -p3 + k, fe)(29) 
E(p) = J dkV(p,k,k)S(k) (30) 

The initial condition is V(p\,p2,Ps) = Vo(pi, P2,Pz) at s = 0. Eq. (29) has a unique 
solution, which satisfies 

V{pi,P2,P3) = V (p 1 ,p 2 ,P3) (31) 



+ 



J dk V (pi,k+p 2 -P3,k) G(k)G(k+p 2 -p 3 ) V(k,p 2 ,p 3 ), 



as can be verified by differentiating with respect to s. Iteration of (31) gives the RPA 
series for V. The selfenergy becomes 

t(p) = J dkV (p,k,k)S{k) 

+ J dkV (p,k,k)G(k) 2 J dlV(k,l,l)S(l). (32) 

The last integral over I gives t(k). By (20), GtG = G-S,so 

±{p) = j dkV (p,k,k)G(k). (33) 
With the initial condition E = for s = this integrates to 

E(p) = J dk V (p, k, k)G(k), (34) 

which is the Hartree equation for E because G = (Q — E) _1 . 

The diagrammatic form of this argument is shown in Figure 2. 

It is obvious that a similar argument works for the summation of the crossed particle- 
hole ladders and the Fock term for the selfenergy. If we resum particle-particle ladders, 
the Fock term vanishes unless we pose a small symmetry-breaking initial condition for 
the selfenergy. In the next section we show that this reproduces BCS theory. 
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Figure 2: The RPA + Hartee resummation. (a) The flow equations (29) and (30). The 
dashed line indicates the s-derivative in (28). (b) The self-consistency equation for V 
corresponding to Eq. (31). (c) The graphical representation of the proof leading to Eq. 
(33): in the first line, (b) is inserted into the graph for the selfenergy from (a). In the 
second line, the inner tadpole graph is replaced by the derivative of the selfenergy by 
(a). Finally, the graphical relation corresponding to (20) is used. 
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3 Superconductivity 



We now use the flow equations (25) and (26) to derive an RG flow where an arbitrarily 
small symmetry-breaking field is put in as an initial condition, and we show that this 
produces a flow to a symmetry-broken state. In the standard Wilsonian picture for 
fixed points, a symmetry-breaking field moves the starting point in a relevant direction, 
and thus changes the fixed point at which the flow ends. 



3.1 A small initial superconducting gap 

The role of the external field that breaks the symmetry in magnets, and which is sent 
to zero to study spontaneous symmetry breaking, is a small external field coupling to 
Cooper pair fields. In the fermionic action, it appears as a gap parameter A . 

In our functional integral representation of the many-fermion model, the fermionic 
fields ip a (k) and ip a (k) depend on the momentum k = (ui, k), where u) is the fermionic 
Matsubara frequency and k the spatial part of the momentum, and the spin index a = 
±. With the notation / dk F(k) =T^2 W L~ d £ k F(ui, k), the quadratic part of the 
action is 

A = J dk (^(fc)(iw - e(k))V(fc) 

+ Ao(k)^(fe)|^(-fc)-A (k)V(fc)|v(-*))- (35) 

Here ip(k)ip(k) = ip + (k)ip + (k) + V>_ (k)ip- (fc), e(k) = e(k) — fi is the dispersion 

relation, and ip(k)eip(—k) = ^ Q , V> Q (fc)£ QQ 'Vv(— k) with e_| = — e |_ = 1 and 

e ++ = e__ = 0. The first term in Ao is the usual Fermi gas kinetic term. The 
dispersion relation e includes the chemical potential fi. The second contains a gap 
function Ao(k) which acts as an external field coupling to the Cooper pairs. The 
Cooper pair is a singlet, so A (— k) = A (k). We also assume that A depends only 
on k. A nonzero A breaks the charge symmetry. As A — > 0, the charge symmetry is 
restored in the action. The question we take up here is what happens with expectation 
values of a model where we add an attractive interaction of the fermions. Spontaneous 
symmetry breaking means that the limit of ( - 0£'0) is nonvanishing in the limit Ao — > 0. 

To make contact with the formulas derived in the previous sections, we introduce a 
Nambu-like field ^ as the column vector 

*(w,fc) - (<A + (fc),^(-fc)>+(fc),TM-fc)) T (36) 

(with the T indicating the transpose) and rewrite the action as Ao — 5(^1 Q^) where 
the brackets denote the bilinear form (/, g) = Y^j=i I dk fj (k)gj (k), with j indexing 
the four components of and the matrix 



A (k) iw-e(k) 

-A (k) -iw-e(-k) 

-iw + e(k) -A (k) 

iw + e(-k) A (k) 



(37) 
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Note that our conventions differ from standard Nambu notation in that the usual off- 
diagonal terms appear in the block-diagonal. This is because we want to have an 
antisymmetric matrix, corresponding to the general formulas given in Section 2. 1 . We 
have also changed notation slightly, denoting by Q what was called Q in Section 2. 1 . 

We are only interested in dispersion relations e satisfying e(— k) = e(k), in which 
case we can write O in the block form 



eA (k) 
-q(fc) 



q(fc) 
-eAo(k) 



(38) 



with 



q(fc) = iuja 3 - e(k), 



e as above, and 03 = diag {1, —1}. 

The fields in position space are defined as 



dk e 



\Vjk-x 



(39) 



(40) 



with x = (r, x) and k ■ x = —lot + k • x. We take (v\, . . . , U4) = (—1, 1, 1,-1) to 
get the convention that -tjj and ip Fourier transform with opposite signs in the exponent. 
This is necessary for the kinetic term to be of the form in (35) because the bilinear 
form (-, •) does not involve a complex conjugate. 



3.2 The RG flow for the reduced BCS model 

The RG scale e s = eq c~ s is used as an infrared cutoff. That is, we replace q(fc) in 
(38) by q(s, k) = q(fc)x(s, k) -1 , where x(s, k) is a smooth function that becomes 
very small when |q(fc)| < and vanishes for q(fc) = 0, and x( s , k) gets close to 1 for 
|q(fc)| » e s (for the diagonal matrix q(fc), |q(fc)| = \\lu — e(k)|; in general one would 
take the smallest eigenvalue of q(fc)). For the proof given in this Section that the gap 
flow leads to the exact solution of the BCS model we do not need to assume a specific 
form of x, but only that x(s, k) — > 1 for s — > 00. In particular, \ may also depend on 
the frequency 10 (recall that k = (ui, k)). For the specific calculations in Section 4, we 
shall take a frequency-independent cutoff function and then take the limit of a sharp 
cutoff. 

In the RG flow, the gap function becomes a function A(s) of the scale parameter s. 
If the interaction is attractive, A(s) increases with s, that is, as the energy scale e s is 
lowered. 

We now derive the flow equations for the reduced BCS model directly from the 
functional form of the equations (25) and (26). We also show that the third order term 
in (25) is indeed irrelevant in this model in the thermodynamic limit. The essential 
reason for this is that the special structure of the BCS interaction fixes loop variables 
in this third order contribution. 
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3.3 The reduced BCS model 



The Hamiltonian of the usual reduced BCS model is of the form of a free fermion 
Hamiltonian H plus a mean-field interaction 

H = H - g L- d &C, C = L- d Y, /(k)c k;+ c_ k ,_ (41) 

k 

The fermion operators Ck, Q have commutation relations {c k a , cj, a , } — L d 8 a ^ a > <5 k ,k' • 
The function /(k) is the gap symmetry function. We have assumed singlet pairing, 
so /(— k) = /(k). We assume that g > so that the interaction is attractive. Al- 
though the reduced interaction contains one momentum sum less than a more realistic 
short-range interaction, the inverse volume in (41) is kept. This implies that in the 
corresponding functional integral, the frequency dependence of the interaction gets re- 
duced to zero frequency pairs in the thermodynamic limit [20], so that the interaction 
effectively becomes 

S r Bcs - -jjXX (42) 

with fl = /3L d and 

uj,k 

X = ^E/( k )^' k )|^(-^-k). (43) 

Again e = ia 2 - The functional integral with interaction (43) can be solved explicitly 
in the sense that it can be reduced to an integral over a single complex variable. In the 
following, we calculate the RG flow for this model. 

The propagator in the flow equation is taken as G = with Q given in (38), but 
where Ao gets replaced by a scale dependent gap A(s). For the interaction, we make 
the ansatz 

V (*) = ~ (vXX + U -XX + + O ((f!)- 2 ) (44) 

where ft — (3L d . The XX and XX terms correspond to the non-charge-invariant 
parts of the vertex. Consequently, the initial condition for the vertices v = v(s) and 
w — w(s) is v(0) — go and w(0) = 0. The coefficient of XX is fixed to be the 
complex conjugate of w to have !7(1) invariance of the initial interaction. 

The RG equations for v, w and A(fc) can now be straightforwardly obtained from 
(25) and (26) by calculating the action of the Laplacians Ac- This is done in Appendix 
A. In that appendix, we also show that our ansatz, in particular (44), is complete for the 
reduced BCS model. That is, f2 plays a role analogous to N in a large- N expansion, 
and the normal part of the selfenergy and any quartic term with a structure differing 
from that in (44) gets an additional factor I/O. In particular, the third order term in 
(25) only produces contributions that are of order Q^ 1 compared to the terms in (44). 
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It turns out (not surprisingly) that the k-dependence of A must be A(k) = f(k)8 
with a scale-dependent S = 8(s) and the function / given in the reduced BCS interac- 
tion. The initial condition for 5 is 8(0) = So with a nonzero So- We take <5 > to fix 
the phase of the gap. Then w is real and the flow runs in the three-dimensional space 
(v, w, 8). The flow equations are 

v = B(v 2 + w 2 ) + 2Avw (45) 
w = A(v 2 + w 2 ) + 2Bvw (46) 
8 = C{v + w) (47) 



with 



A = -^E/ dk /( k ) 2a G( w ' k ) 2 («) 
B = iE|dk|/(k)| 2 ( w 2 + e (k) 2 )6 G (c,k) 2 (49) 



C = dk/(k)a s (c,k) 
and (denoting \Q\ 2 {lj, k) = lo 2 + e(k) 2 ) 



(50) 



a ° = IQ| 2 + X (Ax) 2 (51) 

bG = \Q\ 2 +Vx) 2 (52) 



2xxM 2 

'[|Q| 2 + (Ax) 2 



The term A is the particle-particle bubble with anomalous propagators, B is the particle- 
particle bubble with normal propagators, and C is the loop integral corresponding to 
the Fock term with an anomalous single-scale propagator. 
The equations for g = v + w and — v — w decouple: 

g = (A + B)g 2 (54) 
9 = (B-A)9 2 (55) 

so their solution is 

«"- i-(/1Vo) ,W ° hfW <56) 

Because the initial interaction is charge invariant, g(0) = 9(0) = g - Only the combi- 
nation g = v + w appears in the equation for the gap. By (56), the solutions for v and 
w obey the equations 

v = go + go(Bv + Aw) w — go(Av + Bw). (57) 
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We can now show that the solution 5 of the flow for the gap satisfies the BCS gap 
equation in the limit s — > oo, where the RG scale eoe~ s vanishes. Taking a derivative 
of ac with respect to the scale parameter s we get the equation corresponding to (20), 



a s = a d -A (\Q\ 2 b% - a 2 G ) . (58) 



Thus 



and 



C = \ E / dk /( k ) 2 «c( fc ) - A ( A + B ) 



(59) 



(60) 



A(l + g(A + B))=g^ f dkf(k) 2 a G (k). 
Since the solution for g obeys 

30 = l + d+B)g> m 
we see that the equation for the gap becomes 

^ = 9o^J2 J dk f^G^) (62) 



which integrates to 

P^fJ uj 2 + e(k) 2 + (5(s)/(k)x s (k)) 2 



Here we have made explicit the s-dependence in the notation. In the limit s — > oo, 
X — > 1 and we obtain 

*-*=4E/ dk fl k > ^ + e(k)'%/(k))' (64) 
For <5o — > this becomes the BCS gap equation for 5: 

S - 80 /dk /W ^,^ - ==, t»nh + (05) 

The limiting procedure <5o — > picks the positive solution of this equation. The choice 
of cutoff function x did not play any role in this argument; this is an example for 
universality. 



3.4 The graphical representation of the argument 

In Figures 3 and 4, we show the vertices and flow equations and exhibit the proof 
given in the previous subsection that they resum ladders and anomalous Fock diagrams 
leading to the BCS equation in graphical form. The graphical argument is closely 
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Figure 3: (a) The vertices and propagators in the superconducting gap flow. (b),(c) RG 
flow equations for the normal vertex v (Eq. 45) and anomalous vertex w (Eq. 46), 
respectively. The dashed line indicates that the derivatives of A and B with respect to s 
appear in the flow equations, (d) RG equation for the off-diagonal selfenergy (Eq. 47). 



similar to that depicted in Figure 2 except that there are now two types of interaction 
vertices, and that the relation between S and G involves a matrix, 



G = S - G \ y I G. (66) 




Graphically, A corresponds to a two-legged object with two outgoing arrows and A 
to one with two incoming arrows. For this reason, there are two A terms in the next- 
to-last line of figure 4. One can now write down the equations corresponding to these 
graphs directly from the 1PI graph rules discussed in [4]. 
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Figure 4: The graphical proof that the RG solution satisfies the BCS gap equation. 
(a),(b) self-consistent equations for the normal vertex v and anomalous vertex w, re- 
spectively, corresponding to Eqns. (57). (c) RG equation (47) for the off-diagonal 
selfenergy when (a) and (b) are inserted as solutions for v and to. This yields the 
derivative of the usual BCS gap equation. 
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4 The gap flow 



It is very instructive to look at the solution to the gap equation explicitly because one 
sees how the functions develop when s varies. From this one can read off general 
qualitative properties which we discuss below. 



4.1 Specific choice of cutoffs and the resulting flow equations 



For the explicit representation of the resulting integrals it is convenient to take the 
cutoff function independent of u>. Then the Matsubara sums can be done explicitly. 

A possible choice is x( s >k) = 1 — / 7 (^y l)f-y(— with some 7 ^> 1 and 

f 1 {x) = (1 + c ix )^ 1 . We also take the limit of a sharp cutoff, 7 — > 00, where x 
becomes a step function (for a careful discussion of how this limit is taken, see e.g. 
[11]). Then x 2 — X, the term A% in the denominators in A, B and C simply becomes 
A, and the single scale propagator becomes 



The flow equations then read 

S = gS If 



lo 2 + e(k) 2 + A 2 ' 







9 = 9 



2/e 



= tanh^Vl+W I 
■_ ; j_ s 2 f 2 2 



e = 



dsY 2£ 2 



>e 3 



Here£(k) = ^(k) 2 + S 2 f(k) 2 , 



F(E) = 



1 

2E 



dE \ 2E 2 



( ■ )e s denotes the average on a single scale 

d d k 

<T = ± 



(2tt) c 



S(e(k)-ae s )A(k), 



and (■)>«, denotes the average above scale e s , 



k:|e(k)|>e s 



d d k 

{2ir) d 



A(k). 



(67) 

(68) 

(69) 
(70) 

(71) 
(72) 

(73) 



For an s-wave gap, / = 1 and we can therefore rewrite the equations using the density 
of states N(E). Let 

N(E) = N(E) + N(-E) (74) 
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then 



S = gS N(e s ) — - s tanh £ Je 2 + S 2 (75) 



9 = 9 2 



^ J™ dE N(E) (± tanh ^£ + 25 2 F(£)^ (76) 
d f°° - 1 /3 

— y dEN(E) — tanh^£ (77) 



2£ 2 



where now £ = V-B 2 + 8 2 . 

4.2 Discussion of the solution 

With the sharp cutoff, (65) only gets a restriction \E\ > e s in the integration for finite 
s, so that 

6(s) (l - .go jf° dE N(E)^ tanh - <5 (78) 

From this and (77), we obtain the solution 

at \ - So S(s) S 

e{s) -1-go J~ dE N{E)& tanh f f = ^ ^ ^% ^ 

which is the exact solution for the tangential vertex. It diverges in the limit #o — ► 0; this 
is the remnant of the 1 /q 2 singularity of the Goldstone boson in this reduced model. 
The solution for g is 

5(S) ~ 0(a)-i -2f™dEN(E) S 2 F(£) m 

which for s — > oo converges to the exact solution for the radial vertex. It stays finite in 
the limits s — > oo and A — ► 0. This reflects the fact that amplitude fluctuations of the 
order parameter are gapped. 

The solution of the equations is shown in Figures 5 and 6. Characteristically, the 
gap remains small of order Ao in the flow down to a scale e* . Then it rapidly starts to 
grow and it finally saturates at a value that is the closer to the BCS solution the smaller 
Ao was. For small Ao, the scale e* is only a little bit larger than the scale e c s where 
the flow with A = diverges. For A = 0, the flow of the gap is identical to that 
of the superconducting susceptibility studied in [1, 2, 3, 4]. For A ^ 0, the coupling 
constants become large, but they do not diverge in the flow. 

The qualitative features of the flow for the gap are easily understood from the equa- 
tion (75) itself. The right hand side is positive, so the gap grows with s. On the other 
hand, e s decreases with s, so there is a crossing scale where A(s) = e s . When A s 
gets much larger than e s , the ratio , £s in (75) becomes proportional to e s so the 

ytej+8(s) 2 

scale derivative of the gap vanishes and the gap eventually saturates. This is also true 
for gap functions / that have nodes: in this case, there is always a region of momentum 
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Figure 5: Left plot: Gap amplitude S(s) versus cutoff e s from the RG at zero tem- 
perature, density of states Nq = 1, and initial interaction vo/ e o = 0.3 in units of the 
bandwidth eo. For So = 0, the Cooper instability is at e°/eo = 0.0357. The heavy line 
is for initial gap So /eo = 2.4 x 10~ 4 while the thin solid line is for Sq/eq = 6 x 1CP 8 . 
The dashed line (right on top of the thin solid line) is the asymptotic cutoff-dependence 
S(e s ) = 2^Jt c s {e c s — e s ) for e s < e c s and <5o ^ 0. Right plot: Flow of the normal vertex 
v s (upper curves) and anomalous vertex w s (lower curves) for the two different initial 
gaps. For e s -> 0, w s -> -v s . 



space where 5(s) < e s but this region shrinks to the nodal points for s — > oo and its 
contribution to the right hand side of the flow equation vanishes for s — ► oo. 

The gap flow thus exhibits three regions of clearly distinct behavior, which also 
exist in non-reduced models. Initially, the gap is much smaller than the energy scale, 
so it affects the flow very little and the flow stays close to that of the susceptibility. 
Indeed, as long as the gap is small compared to the energy scale, an expansion of 
the propagators in the gap parameter is convergent because the energy scale e s is an 
infrared cutoff on the propagators. 

In the final stages of the flow, the absolute value of the gap is much larger than the 
energy scale. In a non-mean field model, the gap may still have phase fluctuations. If 
no vortices are present, a gauge transformation can be used to move the phase depen- 
dence from the gap to the hopping term of the effective action, which is much smaller 
than the absolute value of the gap and can therefore be expanded in a convergent series. 
Thus now the gap term is the dominant term in the denominator and it determines the 
long-distance behavior of the fermion propagator. 

The intermediate regime is that where A(s) and e s cross, i.e. the gap is comparable 
to the RG energy scale. In this regime one expects phase fluctuations to be of great 
importance for the detailed behavior of the fermion propagator. 

4.3 Caveats 

As discussed, the general properties of the gap equation already imply that the gap sat- 
urates at low scales. However, they do not imply that the gap saturates at the correct 
value. This happens only because G was used instead of S (the effect of that replace- 
ment on the vertex is even more drastic; see below). Here we show how things go 
wrong when this replacement is not made, i.e. we now consider the truncation to the 
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10" 7 1X»"f 10" 3 10" 1 10" 8 ICft 10" 4 10" 2 10° 

Figure 6: Left plot: Final gap <5fi na i (s — * oo) versus initial gap So- the dashed line 
is the BCS value 5 = 2e exp(-l/iV wo) for v = 0.3 and N Q = 1. Right plot: 
Log-Log-plot of the final coupling strength Ufi na j versus initial gap Sq. The dashed line 
is the BCS result v = Vq6/(28q) for the vertex in presence of a symmetry-breaking 
external field So. The deviations at very small So are numerical errors when the flow of 
the vertices becomes very steep. 

second order flow that is obtained by discarding the 1PI six-point function from the 
flow equations. The only change in the flow equation is that the G gets replaced by 
S. The rest of the structure remains identical. A sample result of such a flow is shown 
in Figure 7. Although the gap grows qualitatively correctly in the flow and saturates 
at a value near the gap given by the BCS equation, it is considerably smaller than the 
correct solution to the gap equation. Moreover, the two-particle vertex diverges at a 
nonzero scale even at positive Ao, in contradiction to the exact solution (79). The 
reason for this behavior is that putting S on the propagator lines instead of G neglects 
contributions to the gap from higher scales. Because the gap term regularizes the prop- 
agator, the speed of the flow of the four-point vertex depends on the gap. Because the 
gap is too small, the vertex grows too fast and has an artificial divergence at a nonzero 
scale. This behavior can be understood in that the Ward identity from global U(l) 
charge symmetry is violated in this truncation. As mentioned, the replacement of S by 
G leads to a rearrangement of the hierarchy that involves all orders; the correct BCS 
flow can be recovered only in this way. 

Finally let us mention that if only the w-vertex is kept, the final gap does saturate at 
small scales, but at a value that diverges in the limit Ao — * 0, so that this approximation 
only reproduces the flow of the superconducting susceptibility in the limit Ao — > 0. 

4.4 The Ward identity 

The Ward identities corresponding to £7(1) gauge transformations of the many-fermion 
system are of central importance for the system, especially in its transition to the 
symmetry-broken phase. The existence and essential properties of the Goldstone bosons 
arising from the spontaneous breaking of the symmetry are directly connected to the 
Ward identities. 

Satisfying the Ward identities is difficult in any approximate scheme. In the flow 
equation approach, finding truncations that lead to the correct Ward identities at the 
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Figure 7: Results in the unmodified 1PI scheme. Left plot: Final gap 5a na }, normalized 
w.r.t. the BCS value, versus initial gap So. The dashed line is the BCS value S = 
2cq exp(— 1/NqVq) for vq — 0.3 and No = 1, Right plot: Log-Log-plot of the final 
coupling strength ugnai (solid line) and the anomalous vertex u>fi na i x — 1 (dashed- 
dotted line) versus initial gap Sq. The dashed line is the BCS result v = vqS/ (2Sq) for 
the vertex in presence of a symmetry-breaking external field So- For values of Sq/^o 
smaller than about 10~ 2 , the integration diverges at a nonzero scale e s . 



end of the flow is a nontrivial problem, which was addressed in [17]. A general inves- 
tigation of the role of Ward identities in functional RG flows is done in [19]. 

In the reduced BCS model, the Ward identity is reduced to that corresponding to 
global gauge transformations. Although this is only a one-parameter subgroup, its 
preservation is nontrivial in approximate schemes: the symmetry breaking relates the 
value of the order parameter with a divergence of the two-particle vertex at scale zero. 
We show explicitly below that this Ward identity is satisfied in the full propagator 
scheme. It is violated in other truncations, in particular in the standard formulation of 
the 1PI flow, where S is not replaced by G. This is the deeper reason for the above- 
discussed failure of the unmodified 1PI scheme in the symmetry-broken phase. 

4.4.1 General structure of Ward identities 

Ward identities can be derived under very general conditions and for general sym- 
metries. We assume that we have a linear invertible transformation U of the fields, 
— >• vp' = which leaves the initial interaction potential Vb and the integration 
measure invariant, 

Vb(U*) = VbW andX>(U*) = (81) 

In general, such a symmetry transformation can mix the fields at different space-time 
points and act nontrivially on all internal indices, i.e. written in detail, it reads 

tf'(X') = J dX U(X',X)V(X) (82) 

A change of integration variable in (1) from ^ to implies 

W{U T QV, V T H) = W{Q, H). (83) 
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Here we have again made the dependence of W on the quadratic part Q of the action in 
(1) explicit. The Ward identity for W can be generated directly from this equation. It 
is already clear from this equation that the Ward identity has a form very similar to that 
of the RG flow equation: the symmetry transformation also changes only the quadratic 
part Q of the action. An infinitesimal change can be rewritten in terms of a Laplacian 
with respect to the source terms in the same way as was used in [4] to derive the flow 
equation. However, now the scale derivative of Q is replaced by a commutator of Q 
with the symmetry transformation and an additional dilatation operator appears, as will 
be explained now in more detail. 

So let us assume further that the group of symmetries U is a Lie group and consider 
an infinitesimal transformation U = 1 + A + 0(A 2 ). To linear order in A we have 

- (AH, ^L) W + c w ^A ATQ+QA e~ w = 0. (84) 

For A T = —A, the operator in the Laplacian is a commutator. In general (and in our 
specific setting) one has commutators and anticommutators; see below. Obviously, the 
form of this equation is similar to that of the flow equation (2) for W. The reason for 
this similarity is that both the scale change and the symmetry transformation change 
only the quadratic term of the action, not the interaction potential Vq. By the standard 
relations between W and its Legendre transform T, we obtain the Ward identity for T 

( A § ' *) - \ Tr A] (0) ') + ( * ' [0 ' A] * } = °" (85) 

In the next subsection, we specialize to the identities we need for our specific case of a 
spontaneously broken U(l) symmetry. 

4.4.2 The U(l) Ward identities for the superconducting state 

To make the effect of the symmetry-breaking terms explicit in the subsequent discus- 
sion, we switch back from the Nambu-type field to the fields ip and jp. Thus 

Z(fj,r],Q,A) = J VipVip e -(^>Q^)-i(^^)+i(^<^) c - v o(^-,4>)+(v,-4>)+(4>,v) 

(86) 

A general gap term is given by 

with an antisymmetric function A of £ = (x, a) and The gauge transformation is 

Vv(x) = e ia ^ ^(s), Mx) = c~ ia ^ (88) 

As before we get 

Z(fj, 77, Q, A ) = Z(e ia f], e^r?, C - ia Qe ia , e ia A e ia ). (89) 
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Taking an infinitesimal transformation, we get for W — 



-InZ 



W + c 



— c 



w - 



(^,i{a,Ao}^) + (^,i{a,A^}A) 



-ii' 



-w 



(90) 



By Legendre transformation, this implies the following Ward identity for the 1PI gen- 
erating function T(<f>, 4>) 



= 



+ 



M,f)-(ia^,f) 



(0, i[a, Q]<t>) - - [(<f>, i{a, A o }0) + {4>, \{a, A }4>)] 
i{a,A } / V J 



1 

2 Tr 



2 

i{a, A^} 
-i[a,Q] T 



(91) 



Again there is some similarity with the flow equation, but now also off-diagonal terms 
appear, because the gauge transformation affects all quadratic terms in the action. Here 



G~ 



A (Q - E) 

-(Q-£) T -A 



(92) 



The usual Ward identities are now generated by expanding in the fields. In second 
order in <f> and <j), only the four-point vertex V = contributes and we obtain a 
relation between the normal and anomalous selfenergy and the corresponding parts of 
the vertex. We shall need here only the relation for the reduced BCS model and for a 
gauge transformation with constant a. In that case [a, Q] = and {a, A } = 2aA . 
Moreover, for the singlet pairing considered here, A((x,a),(x' ,a')) = Ecr ta 'A(x—x') 
where A has Fourier transform A(k) = f(k)S, and in the thermodynamic limit S = 0. 
Comparing the (fxf) term gives the relation 



5 - S = S (v 



dk /(k) 5 



uj 2 + e(k) 2 + ((5/(k)) 2 



Together with the gap equation (64), this implies 

S_ 

T 



v — w 



9o 



(93) 



(94) 



which is the exact solution for the tangential vertex v — w. Thus, given A, the vertex 
is determined by the Ward identity, and the exact result (79) for the vertex is recovered 
only if the gap satisfies the correct gap equation. 



5 The gap flow for general models 

The reduced BCS model is a mean-field model for which the resummation of the 
particle-particle ladders and the Fock equation for the anomalous selfenergy are ex- 
act. We have seen that these features also come out naturally in our gap flow. If we had 
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not introduced the initial gap, we would have found a flow to strong coupling with a 
divergence of the coupling function at a scale just below the one where the gap starts to 
grow in the flow. In this way our calculation also provides an example how a fermionic 
flow to strong coupling is converted to a regular flow to a symmetry-broken state. In 
our example, this state corresponds to the superconducting fixed point of the RG. 

Our method also applies to more realistic short-range models. In these cases, the 
interaction vertex includes more terms than in the reduced BCS model. Then the 
particle-hole diagrams will give nonzero contributions in the flow. Moreover, there 
are anomalous "3:1" vertices with three particles going in and one coming out, as well 
as others that vanish in the thermodynamic limit of the reduced BCS model. With 
these straightforward generalizations, the method can be applied to any model that has 
superconducting correlations. 

Because the gap flow is a smooth modification of the susceptibility flow it au- 
tomatically generates the correct gap symmetry in the course of the flow - wrong gap 
symmetries remain suppressed in the gap flow just as they are in the susceptibility flow. 
One can also study flows where the shape function / of the gap is not fixed but itself 
determined during the flow. Let us give a typical example. In the repulsive Hubbard 
model on the two-dimensional square lattice, the initial interaction does not contain 
any attractive pairing channels. However in the course of the RG flow, particle-hole 
corrections generate an attractive d x 2_ y 2-wave component. Stopping the RG flow a 
certain scale and then resorting to a mean-field solution of the theory at lower scales 
can lead to ambiguities regarding the correct scale for joining these two treatments, and 
all interaction terms involving ipki>k+q with small q ^ have to be discarded. The gap 
flow finds the order parameter without these difficulties. 

It is also possible to include more than one order parameter in the initial conditions, 
hence to study the competition of different ordering tendencies, such as antiferromag- 
netism and superconductivity. In summary, we believe that the method has the potential 
to become a useful and convenient calculational tool for systems with competing inter- 
actions. 

An alternative way of viewing the initial scale A of the gap flow is as follows. 
The initial interaction may contain many other terms in addition to the interaction that, 
either directly or by higher order effects, generates superconductivity. Since the in- 
teraction in a real system is determined only up to some small error terms, one can 
imagine that some very small interactions are present that have a nonzero component 
of a pairing interaction of the form g\XX, where X and X are as in the reduced BCS 
model, with some gap symmetry function /. If we perform a Hubbard-Stratonovich 
transformation on these small terms, the partition function becomes 



Now the HS field h plays the role of the initial gap Ao, and the fermionic integral at 
fixed h can be done by the RG flow. The only difference is that instead of taking the 
limit h — > 0, we now have to integrate over h. The smaller g\, the more concentrated 
the integral will be around h = 0. The point h = itself leads to a divergent flow, 
but in the integral over h, this is inessential because at any other value of \h\, the flow 




(95) 
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converges to a finite result, which gets closer to the BCS gap if \h\ is made smaller (and 
if the gap symmetry is the correct one). In other words, the exceptional point h = is 
really irrelevant by a measure zero argument. 

It is important to note that in this approach, there is still the overall integration over 
the phase ip of h which is part of the h integral. It leads to a vanishing expectation 
value of X even in the symmetry broken phase. This is easily understood: picking 
a particular ip corresponds to selecting a particular pure phase. Taking the average 
over ip produces a mixed phase in which noninvariant expectation values vanish but 
off-diagonal long range order persists. 



A Derivation of the flow equations 

The full propagator is 



lu 2 + e(k) 2 + (xA) 2 V Q ^ A X 
The single-scale propagator is 

<? = * ( -2e A xqq B \ 

^ + e(k) 2 + ( X A) 2 V ~B 2qqeAx J { ' 

where B = q(|A 2 |x 2 - qq) = B T . 

A.l The main terms 

Consider the term quadratic in V on the RHS of (25). The square of the Laplacian 
removes four fields from the degree 8 polynomial in the fields. We can group the terms 
in two categories, which we now discuss at the example of one term in the product 
y (1) y (2) , namely 

'A^X^X^X^XW] (98) 

Here we abbreviated A = A[v 2) and X = X^\ The first type of terms is the 
one where two of the four factors remain untouched by derivatives. There are four such 
terms. Calling the sum of all other terms 1Z, we have 



ww XX 



A 2 (x«x( 2 >) ] +11 (99) 



1Z gets contributions from the terms of the second type where derivatives act on at 
least three factors. All terms contributing to 1Z vanish in the limit L — > oo because 
a loop integration gets fixed by an external momentum (this is specific to the reduced 
BCS interaction and tracing this effect gives a diagrammatic argument that the reduced 
BCS model is exactly solvable). We illustrate it at one example in the next subsection. 
Dropping these remainder terms 1Z in all contributions, calling 

B = ^2IA 2 G (1 (1) 1 (2) )] = ^21A 2 G (X {1) X (2) )] (100) 
A = ±2[A 2 G {XWXW)], (101) 
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and using that 

i2[A 2 G (X«X( 2 ))] = A (102) 

we get the equations 

v = B(v 2 + ww) + (Avw + Avw) 
w = Av 2 + Aw 2 + 2Bvw 
w = Av 2 + Aw 2 + 2Bvw 

If we assume A > 0, w = w is real and these equations simplify to 

v = B(v 2 + w 2 ) + 2Avw 
w = A(v 2 + w 2 ) + 2Bvw. 

The equation for the scale-dependent gap is 

5 = vA s X + wA s X (108) 

(as in the coupling flow, the normal part of the selfenergy vanishes for L — > oo because 
a loop integration gets fixed due to the special structure of the reduced BCS interaction, 
and there remains an overall factor L~ d ). For A > this again simplifies to 

A = gA s X (109) 

The graphical interpretation of the coefficients A, B and C is straightforward: A and 
B denote the particle-particle bubbles of the theory; A is the one with anomalous 
terms and B the one with normal terms. C is the value of the Fock diagram with the 
single-scale propagator S on the line. 

A.2 The remainder terms 

In this section, we first show two prototypical terms which turn out to be one power of 
ft smaller than the main terms, and hence vanishing in the limit — ► oo, relative to 
the main terms. After that, we give the general argument for the reduced BCS model. 
This will also make clear how the terms that combine with S to form G are singled out 
in a natural way in this model, hence we give a specific and complete justification for 
the replacement of S by G in the 1PI equations for this example. 

Recall that the four-point function is of order I/O in the reduced BCS model right 
from the start, so subleading corrections have to be of order 1/fi 2 . The two-point 
function and the selfenergy are of order 1. That is, in the limit Q — > oo, the model 
describes generalized free fields. 

We start with a few preparations. In terms of the field components ip and ip, the 
Laplacian reads 

(6 5 \ ( 8 5 \ ( S S \ 

w ~ qbG W w £aG w + w £aG W (110) 
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(103) 
(104) 
(105) 



(106) 
(107) 



(a) 



(b) 



Figure 8: Examples of contributions that vanish in the thermodynamic limit of the 
reduced BCS model, (a) the lowest order normal selfenergy. (b) A second order term 
with three fermions going in and one coming out. 



and similarly for S. Here we used 

qb G ea G 



G= | ~!? G ]. (!!!) 



The functional derivatives obey 

5 



5ip a (k) 



^ a >(k f ) = S a>a > 5 a (k,k') (112) 



where Sn(k, k') — f38 w ^i £ d ^k,k'- With this convention, one can take the limit L — > 
oo that we are interested in, as well as the limit T — > 0, without any further rescaling 
of the fields. 

The first example is the contribution to the right hand side of the flow equation for 
the selfenergy, eq. (26), from the normal Fock graph shown in Figure 8 (a). This is 
a term that changes the normal part of the selfenergy, i.e. it shifts the Fermi surface 
and changes the quasiparticle weight. Its value is obtained by applying the part of the 
Laplacian to the vertex, 

r =-(^ 6 4) h* x - (n3) 

Straightforward evaluation of the derivatives gives 

k.a 

which is of order f2 _1 as compared to the term in the action, hence vanishes for f2 — > 
oo. 

The second example is the contribution to the right hand side of the flow equation 
(25) for the anomalous part of the vertex where three particles go in and one comes 
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out, shown in Figure 8 (b), with value 



' <5^< 2 > 

Straightforward evaluation of the action of the Laplacian gives 

, 2 



(115) 



(116) 



with 

Y 



ds 



(b G (k)a G (k)) 



(ic- e (k))^^ 2 )(fc)^ 2 )(fc). (117) 



Like X and X, Y is a bilinear in the fermion fields with a similar momentum structure. 
However, the prefactor in (116) contains an additional factor 1/fi as compared to the 
main term. 

In both examples, the extra power 1/51 arises because there is no summation over 
loop momenta - the structure of the graphs and the vertices of the reduced BCS model 
fix the loop momenta in terms of the external ones. Then there remains an overall 1 /f2 
which makes the terms subleading in fl 



A.3 The general argument 

It is now easy to generalize this argument to general graphs. The Cooper pair terms X 
and X can be rewritten as 



k.i 



(118) 
(119) 



k',V 



which makes it clear that quantities of order 1 contain one sum over momenta and 
frequencies per field t/j or f/S, if all momentum constraints are written in terms of the 
delta functions Sq. By definition of 6q, 



(S n (k,k')f = Q6 Q (k,k') 



(120) 



Because the initial vertex function is 0(fi _1 ), a term of orderp in the vertex is of order 
0~ p+£ where I is the number of times two delta functions get paired so as to give a 
square as in (120). To get a contribution of order O -1 to the two-particle vertex, I must 
equal p—1. Such a complete pairing occurs only in ladder graphs. Similarly, the 0(1) 
contribution to the selfenergy can come only from anomalous Fock-type diagrams for 
the selfenergy made by joining two external legs of a ladder graph. 

The argument we just gave does not require that the model is doubly reduced. All 
these graphical arguments apply directly to the reduced BCS model in Hamiltonian 
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form, i.e. with full frequency dependence. It should be noted that the diagrammatic 
approach, i.e. determining for each diagram separately whether its contribution is O(l) 
or subleading, is not a proof in the mathematical sense that this gives the exact solution 
because a graph-by-graph analysis presupposes that the perturbative series is conver- 
gent or at least asymptotic (the functional integral argument in [20] also involves an 
unjustified exchange of the thermodynamic limit with the limit of continuous Euclid- 
ian time). A complete mathematical proof that this gives the solution of the reduced 
BCS model in Hamiltonian form was given in [21]; see also [22]. 
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